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I. INTRODUCTION 

The 0{N) model with the spontaneously bro- 
ken symmetry is a traditional touchstone for new 
theoretical approaches with applications to a vari- 
ety of physical phenomena, such as the chiral phase 
transition in nuclear matter the formation of 

Bose-Einstein condensates in atomic gases 0, 0| , 
or inflationary cosmology scenario However, 
even for such a simple model, conceptual problems 
arise for the dynamic and even thermodynamic 
treatment of the system. To be precise, the present 
discussion concerns non-perturbative approaches 
based on partial resummation schemes of various 
kinds: the ^-derivable approximation of Baym 
or CJT formalism Q as it is referred in the field 
theory, 1/7V expansions constructed on top of the 
CJT formalism 0,^^, gapless approximations (34j 
approximations of mixed type like the vari- 
ational two-particle-point-irreducible (2PPI) con- 
cepts'^ [13 or the Baym-Grinstein approxima- 
tion ^ 

Central objective of this paper is to construct 
self-consistent dynamical equations of motion for 
the mean-field (j){x) and one-particle propagator 
G(x,y) with improved features with respect to 
symmetry properties. The space-time changes of 
both quantities can be formulated as a coupled set 
of mean-field and Schwinger-Dyson (SD) equations 
in their differential form 

G^^Qm = J[<j,,G], (1.1) 
G^^ QG{x,y) = I + Y.[4),G]QG. (1.2) 
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Here Gq ^ is the inverse free propagator, the unity 
I and the 0-opcration are to be understood in 
the corresponding functional sense. The self- 
consistence of the scheme requires both source 
terms, the current J and the self-energy S, to be 
entirely functionals of (f) and G. Baym |7l Il5l|. see 
also [13 , has shown that such a scheme is conserv- 
ing at the expectation value level and at the same 
time thermodynamically consistent if and only if 
both driving terms are derived from a single func- 
tional (f>[(/), G] through functional variations 

J=^ci,[^,G], (1.3) 

()Ct21 

Approaches of such <f>-derivablc type permit trun- 
cations of the ^-functional at any level, this way 
providing a variety of approximations without vi- 
olating the above features. Equipped with the re- 
quired irreducibility features, i.e. that the closed 
diagrams of $ are two-particle irreducible (2PI), 
any such truncation scheme, through the coupled 
set of equations l|l.l|l and (|1.2|l . provides a partial 
resummation of diagrams to infinite order, which 
is void of double counting and with built-in de- 
tailed balance. In its Wigner transformed ver- 
sion it leads to the well known Kadanoff-Baym 
equations which through gradient approxima- 
tion permits the derivation of generalized transport 
equations [l7llllll9l|. 

Despite all the above features, however, <I>- 
derivable approximations in the sense of Eqs Hl.l|l - 
1)1. 4|l violate Ward-Takahashi identities beyond 
the one-point level. In particular, this leads to a 
violation of the Nambu-Goldstone (NG) theorem 
IQQiHOl in case of spontaneously broken symme- 
try (which is referred to as the Hugenholtz-Pines 
theorem in the physics of Bose-Einstein condcn- 
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sation). The standard remedy is to construct a 
IPI effective action F [Cj [(/>], (j>] based on the cho- 
sen ^-derivable approximation, i.e. G$[(?!)] being a 
solution of Eq. Hl.l|) with arbitrary field (j>. The 
functional variation of this effective action with re- 
spect to (j) results in equations for higher order ver- 
tices, e.g. the Bethe-Salpeter (BS) equation. This 
BS-equation (the RPA equation in the case of the 
HF approximation) has to be solved supplemen- 
tary to the SD equation with frozen propagators of 
the SD scheme (HH), cf. Ref. The thus con- 

structed vertex functions, in particular the RPA 
self energy, do indeed respect the corresponding 
Ward-Takahashi identities and hence the NG theo- 
rem. However, these quantities do not take part in 
the self-consistent scheme in the sense that they do 
not enter back into the source terms in Eq. Hl.l|) . 
In turn, the internal propagators entering in 
arc still not gapless and still suffer from violating 
the NG theorem. In view of non-equilibrium appli- 
cations, where all dynamical quantities are evolved 
in time, this is an unacceptable deficiency of the 
remedy scheme via higher-order vertex functions. 
Moreover, the corresponding BS equation numer- 
ically requires huge computing resources once one 
includes higher order diagrams in the $ functional. 
Generally this remedy scheme, constructed "on- 
top" of a self-consistent (f>-derivablc SD approxi- 
mation, is of a "gap less" type in terms of classi- 
fication of Ref. [ij. This gapless approximation 
respects the NG theorem, yet violates conservation 
laws and thermodynamic consistency 

In this respect, the particular statement in Ref. 
|2ll| and some subsequent publications that "the 
NG theorem is always satisfied at any finite TV 
in the Hartree-Fock approximation" of the $- 
functional (or CJT) formalism, is certainly mis- 
leading. This statement is based on the identifi- 
cation of the pion and sigma masses with the cur- 
vatures of the effective potential at its minimum. 
These curvature masses are actually precisely those 
given by the RPA scheme. They are again con- 
structed "on top" of the self-consistent scheme as 
they nowhere re-enter in the driving terms of the 
coupled meand-field-SD scheme (ll.l|) - (|1.2|) . What 
however matters in self-consistent dynamical treat- 
ments, both for transport and thermodynamics, 
are the pole masses of the self-consistent Green 
functions. This fact was realized long a go i n the 
theory of Bose-Einstein condensation [sl lllj. As 
a consequence, kinetic descriptions of condensed 
systems were based on the gapless approximations 
[3II11II rather than the ^-functional, accepting the 
lack of conservation laws as a minor defect com- 
pared to the else poorly reproduced Goldstone bo- 
son. 



The above statement can also formulated in for- 
mal terms. Strictly speaking, the NG theorem [2^ 
claims that the inverse propagator of the pion van- 
ishes as as 0. If the Green function is 
calculated as the second variation of the effective 
action over fields (like in the "gapless scheme"), 
then the zero curvature of the effective potential 
at its minimum in the "pion direction" is equiv- 
alent to the the NG theorem. However, this is 
precisely the problem that in the ^-functional (or 
CJT) formalism the Green function is associated 
with stationary point of the effective action with 
respect to the variation over the Green function G 
itself rather than the fields. In this case, the fact 
that the curvature of the effective potential at its 
minimum is zero is irrelevant to the NG theorem. 
This is why the recovery of the "gapless scheme" 
was the main idea behind the restoration of the 
NG theorem, reported in Ref. [23|. In fact, similar 
arguments on the validity of the NG theorem have 
been already put forward in Ref. . 

From this point of view, the two-particle-point- 
irreducible |35l| (2PPI) approach [13 still violates 
the NG theorem. In leading (zero) order in the 
1/N expansion the model is indeed free of the 
above-discussed problems. However, with N = 2 
for Bose-Einstein condensation and iV = 4 for 
the chiral phase transition in nuclear matter this 
zero order scheme is far from being adequate. Al- 
ready terms of next-to-leading order in I /N in the 
self-energies violate the NG theorem in the above 
sense, i.e. due to the finite pole mass of the pion, 
contrary to the statement of Ref. This is 

clearly seen in the thermodynamic limit of this for- 
malism [2^ . 

First attempts to find a compromise between 
derivable and gapless schemes were undertaken by 
Baym and Grinstein In fact, their modified 

Hartree-Fock approximation is a simplified gap- 
less scheme, which still respects the NG theorem. 
Therefore it was accepted as a main tool for stud- 
ies of disoriented chiral condensates [23. 125I |2^ . 
Baym and Grinstein £^lso found severe prob- 
lems with the renormalization of partial resumma- 
tion schemes. Great progress in the proper renor- 
malization of such schemes was recently achieved 
in Refs ^ [HIH IH [s^ . However, we will not 
touch the question of renormalization in this ar- 
ticle, since it leads to extra complications which 
deserve a special discussion. 

In this paper we discuss the standard exam- 
ple of a spontaneously broken 0{N) model in 
derivable Hartree-Fock approximation |3^. The 
considerations are confined to the thermodynamic 
equilibrium, since the discussed problems reveal 
themselves already at that level. Our goal is to 
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construct a modified Hartree-Fock approximation, 
which is stiU ^-derivable, and thus automatically 
provides conservation laws and thermodynamically 
consistence at the level of the corresponding mean- 
field-Schwinger-Dyson scheme and at the 

same time respects the NG theorem. The present 
treatment is based on a naive renormalization, 
where all divergent terms are simply omitted. 



II. CONVENTIONAL HARTREE-FOCK 
APPROXIMATION 

We consider the 0(A^)-model Lagrangian 



H -6. 



(2.1) 



where (f> ~ (0i, 02, (/>Ar) is an iV-component 
scalar field, (f)^ = (f)a4>a with summations over a 
implied. For H = this Lagrangian is invariant 
under 0{N) rotations of the fields, li H = 
and < 0, the symmetry of the ground state 
is spontaneously broken down to 0{N — 1), with 
A^ — 1 Goldstone bosons (pions). The external field 
H ■ 4> = Ha(t>a is a term which explicitly breaks the 
0{N) symmetry. It is introduced to give the phys- 
ical value of 140 MeV to the pion mass. 

The conventional Hartree-Fock approximation 
to the 0(N) model is defined by the $ 
functional I 



$HF 



4^ 



(2.2) 



where the crosses denote the classical field cj). 
Within the ^-derivable scheme the r.h.s. of the 
equations of motion for the classical field (J) and 
the Green function (self-energy E) follow from the 
functional variation of $hf with respect to the 
classical field (f> and Green function G, respectively 

Ucp + m (p ~J = — J — 



G 



HF 



D- 



5$ 



HF 



(2.3) 



■^HF 



SG 
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where D is the free propagator. In the sponta- 
neously broken-symmetry phase, i.e. for H ~ 
and < 0, the solutions of equation set H2.3|l vi- 
olate the NG theorem, as it has been demonstrated 
in numerous papers (sec e.g. Refs Elllllsl). A 
detailed analysis of these equations in a notation 
similar to ours has been given in Ref. |3ll |. 

It is possible to remedy the defect of the 
Hartree-Fock approximation by correcting the self- 
energy by standard random phase approximation 
(RPA) technique as follows [20llp| 



-'gaplcss 



Q 



AS 



RPA 



(2.4) 



where 



AS 



RPA 




(2.5) 



This AS'^^^ is constructed in terms of Green func- 
tions Ghf obtained from solving the self-consistent 
Hartrec-Fock e quat ions (|2.3|) and in fact represents 
the RPA series [S^. Actually, this correction con- 
verts the scheme into that of gapless rather than of 
<I>-derivable type. This results from the loss of the 
self-consistency, since the RPA self-energies do not 
enter the SD-equation (|2.3() . The most important 
ingredient of this scheme is the pion self-energy 



AS 



RPA 



(2.6) 



{iG-jriGa 



2A 

N 



(iG^iG,) 



which provides the NG theorem. Here 



(iGTriGo 



= / d\iGAk + q/2)iG,{k-q/2i: 
J0 



(2.7) 



where the Matsubara integration 



n— — oo 



(27r)3 



(2.8) 



/(27rinT,q), 



is implied with T being a finite temperature. At 
= 0, by means of the Hartree-Fock equations of 
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motion AY.^P^{k = 0) is reduced to |32 



AE«^^(fc = 0) 



2A 

"77 



{Qn 



(2.9) 



where 



d^kGa{k) 



(2.10) 



is the tadpole diagram of the Green function 



Ga(fc) 



1 



Mr 



(2.11) 



which is the solution of the Hartree-Fock equa- 
tions of motion (|2.3|l with Ala being the effective 
mass of the particle. Upon substitution to H2.4|l . 
this AS^P^ provides G;;(g^pj^^^^ (fc 0) = or 

A^^(gapicss) = + ^^fHk = 0) = in the 
broken-symmetry phase. This is precisely the re- 
quirement of the NG theorem. 

However, the recovery of the NG theorem is 
achieved here on the expense of loss of self- 
consistency (i.e. ^-derivable nature) of the ap- 
proximation. It means that the dynamic treatment 
based on scheme (|2.4|l will not possess proper con- 
servation laws and thermodynamic consistency is 
lost. The latter means in particular, that the num- 
bers of particles calculated as a derivative of the 
thermodynamic potential f2 



1. 



;f}{</),G} =/o(<^) + -Tr(lnG-i) 

^iTr(SG) + <i>{0,G} 



(2.12) 



with Io{4') being the free classical action of <j> field, 
over chemical potential |38l|. on the one side, and 
as the Green function G integrated over momen- 
tum, on the other side, will be different. Moreover, 
the spectral function of the sigma meson calculated 
within this approach has the wrong threshold be- 
havior [23|, which however docs not affect the va- 
lidity of the NG theorem. It is not gaplcss because 
of the still nonzero pole mass of the pion Green 
function H2.11|l . Moreover, the chiral phase transi- 
tion still proves to be of the first rather than of the 
second order, as required by the universality class 
of the 0{N) model. 



III. GAPLESS HARTREE-FOCK (GHF) 
APPROXIMATION 



For the exact theory, i.e. when all diagrams of 
the $ functional are taken into account, both the 
gapless and <i>-derivable schemes become identical. 



It is desirable to have them identical also at a cer- 
tain approximation level such as the Hartree-Fock 
approximation. If not completely identical, then 
at least identical in the Goldstone-boson sector, 
which is of prime importance. The latter means 
that corrections to the pion sclf-encrgy H2.6() pri- 
marily need to be incorporated in a <i>-derivable 
way. In fact, we need far less to satisfy the NG the- 
orem, namely only that the pion self-energy van- 
ishes at vanishing four-momentum p. 

For this purpose one can introduce a phe- 
nomenological symmetry-restoring correction to 
the (f> functional of Eq. H2.2|) 



(A^- 1)A 
2N 



(3.1) 



which precisely produces the required correction 
term, namely the RPA term Ai;^P^(fc = 0) to the 
pion self-energy, see Eq. H2.9|l . This A$ can be 
presented in a manifestly 0(A^)-symmetric form 



A$ 



A 

'2N 



(NQabQba - QaaQbb) , (3.2) 



where summations over a and b are implied and 



Qab= [ d^kGabik) 



(3.3) 



is the tadpole diagram of the Green function Gab- 
Thus, the $ functional taken as a functional of the 
classical field </) and the Green functions Gab for 
the gapless Hartree-Fock approximation is defined 
as 

$gHF[0, Gab] = $Hf[0, Gab] + A$[Gafc] (3.4) 

In the broken-symmetry phase with non-vanishing 
classical field (j) — {4'a} and in an arbitrary, 
not necessarily diagonal, representation the Green 
functions. Gab takes the form 



Gab 



-Gc 



0a 06 



G„ 



(3.5) 



in terms of the pion and sigma propagators, cf. 

e.g-ia. 

The advantage of the 3>gHF is that the resulting 
<I>-derivable approximation obeys the NG theorem, 
at the same time keeping all the pleasant features 
of the (^-derivable approach. Moreover, $gHF does 
not change the Hartree-Fock equation for the clas- 
sical field (/) (|2.3|l . as A$ does not depend on 0, 
cf. Eq. (|2.3|) . This is a desirable feature, since 
the classical field equation coincides in both the 
<I>-derivable and the gapless schemes, and there- 
fore it is reasonable to keep this unchanged also in 
the gaplcss Hartree-Fock approximation. Another 
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advantage of A$ is that it does not change the re- 
sults in the 0(iV)-symnietry restored phase, where 
Ml = Ml and hence = Qa and AE^ = 0. 
Indeed, this phase is described quite reasonably 
within conventional Hartrcc~Fock approximation 
and hence requires no modifications. For the fur- 
ther discussion of this modified <i>-derivable ap- 



proximation we use the notion in terms of CJT 
effective potential, see e.g. 0,1^1^, in order to 
comply with numerous previous considerations in 
the literature. 

The manifestly symmetric form of the CJT ef- 
fective potential in the conventional Hartree-Fock 
approximation reads 



14jf(</>, G) = 2"^ 



2 j,2 



4A^ 



H-(j)+l [ d'^k\ndetG~\k) 
2 J0 

Gba{k) - 1 



+ m^) Sab + ^ {(jy^Sab + '2(t)a<f)b) 



+ {QaaQbb + '^QabQba) , 



(3.6) 



e.g., cf. [23|. The phenomenological symmetry- tadpole, cf. \d.3\i . Here and below, summation 
restoring correction, corresponding to A<1> of Eq. over repeated indices a,b,c,... is implied, if it is 

not pointed out otherwise. All the quantities are 
symmetric with respect to permutations of these 
A,,,„ „ „ , — , indices. 



(EH), is 

AV 

where 



2N 



[NQabQba - QaaQbb) , (3.7) 



The equations for Gcd, i-c for the corresponding 
(3.8) tadpoles Qab, and the fields 0c result from varia- 
tions of VgHF = ^HF + AV over Gcd and (j)c, re- 
is the Green function, and Qab is the corresponding spectively, 



G-,\k) = k'Sab + Ml, 



Mia = m'^Scd + - [0'5cd + 20c0d + SQaa5cd + 2(1 - N)Qcd] , 

He = m^ctyc + [(jy'^ct'c + Qaa(l>c + 2Qcd4'd] ■ 



(3.9) 
(3.10) 



These equations arc in an arbitrary non-diagonal representation. Introducing projectors on tt and a states 

1 



N - 



- {Scd - 0c0d/0^) 



n 



cd 



ac9dl 



one arrives at 



Ml n^^Ml, 



N 



Oa(Pb 



iba\ 



Ml nS^A/f^ = + ]^ t^*^' + + - ^) ('^a06/0') Qba] 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



from Eq. H3.9|l . In order to project the mean- field equation (|3.1U|I on the a-direction, we just multiply it 
by (j)c 



jm^ + + Qaa + 2 {cPaCPb/<f') Qba] | = i?c 



(3.15) 
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In the diagonal representation {(per 7^ 0, iJ^ = H 
and = = 0) these equations take the follow- 
ing form 



Mi 



Ml = 



A [302 + (5 _ 2N)Q, + 3(iV - 1)Q,] 
- A [202 + 2(7V-l)(Q,-Q,)] (3.16) 
A [02 + + (VV - 1)0,] (3.17) 



H 



N 



+ 'iQa + {N -l)Q^) 



(.3.18) 



Here we used Q„ — Qaa and Qt^ = Q-mx in terms of 
definition (|3.3() . The tadpoles consist of two parts: 



)(div) 



(div) 



one of them 

d^k 1 1 

(27r)3 (fc2 +A/2)i/2 2' 



(3.19) 



is divergent. Here we adopt a naive renormaliza- 
tion, which consists in omitting this divergent con- 
tribution. Therefore, the tadpoles are defined as 
follows 



Qa = Ql = 



1 



(27r)3 (fc2 +M2)i/2 

1 

'exp[(fc2 + M2)i/2/T] 



(3.20) 



From these equations it is evident that the NG the- 
orem is fulfilled. Indeed, in the broken-symmetry 
phase [H = 0) the square-bracketed term of the 
field equation H3.18|l equals zero, which precisely 
determines the pion mass, cf. Eq. H3.17|l . 



A. Vacuum (T = 0) 

At r = 0, the quantities under investigation are 
"experimentally" known pj: M^{T = 0) = = 
139 MeV, M^{T = 0) = 600 MeV, and the 

pion decay constant 4>q = fn = 93 MeV. These 
known quantities should satisfy Eqs. H3.16|l - H3.18|l 



at T = 0, i.e. with = Qtx = 0. In order to make 
this set of equations consistent, we should put 



(3.21) 



i.e. precisely the same as at the tree level. Then 
equations (|3.17|) and (|3.18|) become identical. Re- 
solving Eqs. (|3.16|) - (|3.18|) with respect to and 
A in terms of "experimental" quantities m.,^, nic 
and /tt, we arrive at 



A = 



N {j 



2f2 



(3.22) 



A 



-m- = -m;-f-/;. (3.23) 
B. Symmetry Restoration Point Ti (at H = 0) 



li H ~ 0, above some temperature the 0{N) 
symmetry is restored. First, we apply the strong 
condition 



Mi = Mi 



0. 



(3.24) 



These equations determine the symmetry restora- 
tion temperature Ti. In this case all three equa- 
tions H3.16() - (|3.18() reduce to the same single one 
with the solution 



T2 



12 



{N 



-f 



(3.25) 



Let us consider the region in the vicinity of Ti, 
assuming |T - Ti | < Ti, < Ti and Ma < Ti. 
Here we use brief notation 

m2 = -m^lT"^ > 0, s = 0/T, Ma = Ma/T, 



^=-™^ + ^(^^ + 2)l.A^(^ + 2. 



r2 _ r2 



Keeping only linear in Ma terms in Eqs (|3.16|) - 
(|3.18|) accordingly to (|A.3|I . we arrive at the fol- 
lowing set of equations 



2s' 
A 



l{N-l) 



4tt 



= R+ — \ s' - — 



N 



1 r„,~. 



Air 



= s 













y An I 



M„ - M^ 
3M„ + {N- 1)M, 
Mia + {N - l)M^ 



(3.26) 
(3.27) 
(3.28) 
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FIG. 1: Schematic behavior of the solution in the vicinity of the Ti (left panel) and T2 (right panel) points. 



The solution to this set with s = is trivial 



s = 0, 

N An 

_ tt T^- 
~ 3 T2 



(3.29) 



and exists only at T > Ti. This is the phase of 
restored symmetry. If s 7^ 0, then = 0. Hence 
we arrive at the solution 



R 



iV2(iV- 1) 
T 2iV+l 



T2 



2N + 1 12 

A 2A^+ 1 

\(« + 2)l!^, ,3,31) 



2A^+ 112 



implying that there are no solutions with > 
and Ma > simultaneously. In fact, this solution 
exists only at T > Ti, since should be posi- 
tive, as required by definition of Q^. However, at 
T > Ti, the solution is unphysical, since be- 
comes negative: < 0. Of course, there exists a 
finite solution with > and > at T < Ti , 
as it is numerically demonstrated in the next sec- 
tion. This behavior of the solution in the vicinity 
of r = Ti is schematically demonstrated in Fig. 
1. It implies that the phase transition at Ti, if it 
ever occurs, is of the first order, since the physical 
quantities reveal jumps at this point. 



C. Partial Symmetry Restoration Point T2 

(at H = 0) 

The transition point T2 is determined by the 
weak condition 



Mi 



0, Alj ^ 0. 



(3.32) 



The symmetry is not completely restored at this 
point, since the meson masses still remain different. 
A calculation very similar to that above gives 



T4 



[N + 2) 



1 



(7V-1) 



(3.33) 



Strictly speaking, this is not a solution for T2, 
since the r.h.s. of H3.33|) still depends on T2 
through M^(T2). However, it numerically occurs 
that Ma{T2) ~ ttIct, which almost removes this de- 
pendence, and hence 



rp2 

^2 



12 



J TT 



(3.34) 



The behavior of the solution in the vicinity of the 
T2 is schematically demonstrated in Fig. 1. Such a 
behavior implies that the phase transition at T2, if 
it ever occurs, is of the second order, since quanti- 
ties reveal no jumps and only their derivatives are 
discontinues at this point. 



IV. RESULTS FOR iV = 4 

The actual structure of the solution to Eqs 
l|??TC|l - (|TTH|) turns out to be even more involved 
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FIG. 2: Meson masses (left panel) and <j} (right panel) as functions of temperature for the — case. Stable 
and metastable branches are presented by solid and dashed lines, respectively. 



than qualitatively discussed above. This is demon- 
strated for the case of exact 0(4) symmetry, i.e. 
ir = 0, in Fig. 2. 

First of all, there are two different solutions to 



the set of Eqs. l|TTB|l - ljTT5|) . To decide which 
branch is more stable, we should compare the cor- 
responding values of the C JT effective potential, cf. 
Eqs H3.6|) and 13.7|l . which in explicit form reads 



+ A [(5 _ 2N)Ql + {N- XfQl + 6(iV - \)Q„Q^\ , (4.1) 

where 

I 



and Qa are given by Eq. (|3.2U|) with Ma being so- 
lutions of equations of motion. Note also that this 
effective potential is nothing else but the negative 
pressure 

P^-y^YLV- (4.3) 

A stable branch should correspond to a mini- 
mum of the effective potential with respect to 
the field variation, i.e. the second derivative 
rf^V^HF(0, G(0))/d0^ should be positive. As it is 
seen from Eqs (|B.3|l - (|B.5|l . all branches with pos- 
itive masses M^- and Af^ correspond to minima. 
This is precisely the case for all branches discussed 
here. This minimum is only very shallow, when 
Af^ = 0. Therefore, we can only choose between 



stable and metastable branches. The most sta- 
ble branch should have minimal absolute value of 
V^HF, i-e. maximal pressure. 



A. Exact 0(4) Symmetry 

The results for stable and metastable branches 
in the case of exact 0(4) symmetry, i.e. H = 0, are 
displayed in Figs. 2 and 3. In fact, it is impossible 
to resolve these branches near their crossing point 
at Tcross — 360 MeV. However, these branches can 
be distinctly resolved for the case of the partially 
violated symmetry, see Figs. 4 and 5, which con- 
firms our present identification. 
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FIG. 3: Pressure (left panel) and pressure difference between stable and metastable branches (right panel) as 
functions of temperature for m,r = case. The pressure for the stable and metastable branches are presented by 
solid and dashed lines, respectively. 



The stable branch exists at all temperatures. It 
starts at T = from the physical vacuum val- 
ues of the masses and classical field and crosses 
the metastable branch at Tcross- In terms of the 
pressure, they are touching rather than crossing 
(see left panel of Fig. 3). Therefore, no transi- 
tion from one branch to another occurs at Tcross- 
In the broken-symmetry phase, the stable branch 
behaves precisely in accordance with our expecta- 
tions - the pion mass equals zero. Then a phase 
transition of the second order occurs at T2 ~ 180 
MeV, at which the field becomes zero. However, 
the TT and a masses still differ beyond this transi- 
tion point. They become equal only after the sec- 
ond phase transition also of second order at Tcross- 
Note that the equal-mass solution above Tcross is 
precisely the same as in the conventional Hartree- 
Fock approximation, since the gapless modification 
term vanishes in this case. 

The Ti point proves to be irrelevant for the 
stable branch. Rather it is the starting point 
for the metastable branch, which in the range of 
Ti < T < Tcross precisely coincides with the solu- 
tion of the conventional Hartree-Fock approxima- 
tion. The existence of this branch above Tcross is 
an artifact of the introduced NG-thcorcm restoring 
term A$. The corresponding field always vanishes 
for this branch. 

Summarizing, in the spontaneously broken- 
symmetry phase the NG theorem is indeed ful- 
filled. Thus, we achieved the goal of the proposed 
modification of the Hartree-Fock approximation, 
however partially, since the phase transition now 
proceeds through a sequence of two second-order 
phase transitions rather than a single one. 



B. Approximate 0(4) Symmetry 

For the explicitly broken symmetry case (finite 
pion mass) the results for both branches are dis- 
played in Figs. 4 and 5. In this case these two 
branches are well resolved in the entire range of 
temperatures, since they do not even touch each 
other, see Fig. 5 (right panel). As it was ex- 
pected, the sequence of two phase transitions is 
transformed here into smooth cross-over transi- 
tion. 



V. CONCLUSION 

A modified ^-derivable Hartree-Fock approxi- 
mation to the 0(A^)-model is proposed, which si- 
multaneously preserves all the desirable features of 
^-derivable scheme (i.e. conservations and ther- 
modynamic consistency) while respecting the NG 
theorem in the broken-symmetry phase. This 
is achieved by adding a correction A(f> to the 
conventional $ functional. With this correction 
term, the chiral phase transition proceeds through 
a sequence of two second-order phase transitions 
rather than a single one. In the first transition 
the field disappears but the meson masses still re- 
main different. In the second transition also the 
masses become equal, and the 0{N) symmetry is 
completely restored. 

The nature of this correction can be understood 
as follows. For the full theory, i.e. when all dia- 
grams in $ functional are taken into account, the 
gapless and ^-derivable schemes are identical and 
both respect the NG theorem. The conventional 
^-derivable Hartree-Fock approximation omits an 
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FIG. 5: Pressure (left panel) and pressure difference between stable and metastable branches (right panel) as 
functions of temperature for m-n = 139 MeV case (stable and metastable branches presented by solid and dashed 
lines, respectively). 



infinite set of diagrams which are necessary to re- 
store its equivalence with the gapless scheme. The 
A$ correction to the Hartree-Fock approximation 
takes into account a part of those omitted diagrams 
(at the level of actual approximation), and thus re- 
stores this equivalence in the pion sector. 

This A$ correction is unambiguously deter- 
mined proceeding from the following requirements: 
(i) it restores the NG theorem in the broken- 



symmetry phase, (ii) it docs not change the 
Hartrec-Fock equation for classical field, since the 
^-derivable and gapless schemes provide the same 
classical-field equation already without any mod- 
ifications, (iii) it does not change results in the 
phase of restored 0{N) symmetry because there 
is no need for it. If two last requirements are re- 
leased, the A$ correction become ambiguous: any 



^^alternative 2iV L^^^ ^<y) 

f {4>y{N - 1) + 0. - Q., - 1) + 0. + Qa/iN - 1) + const)] , (5.1) 
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where / is an arbitrary function satisfying 
df{x,y)/dx = at a; = 0, is as good as (|3.1(l . 
The arguments of the / function are taken to be 
0{N) invariants. For the first argument it was 
demonstrated in Eqs. (|3.1|1 and (|3.2|) . whereas for 
the second argument it foUows from the identity 
Qaa = Qa + {N — l)QTr. Notc that our actual 
choice (Ell) is as close to the RPA result [H IH 
as possible. 

The treatment of this paper was based on a naive 
rcnormalization scheme, when all divergent parts 
of diagrams were simply omitted. Presently, con- 
sistent rcnormalization schemes of ^-derivable ap- 
proximations are available IHIll H^l ■ Our 
first experience of dealing with them indicates that 
they give rise to new (yet not always desirable) 
features of <i>-derivable approximations to theories 
with spontaneously broken symmetry. These new 
features deserve special studies. To avoid mixing of 
effects of the NG-theorem restoring term A$ with 
those of the renormalization, we defer a discussion 
of the latter to a subsequent paper. 

The introduced gapless Hartree-Fock approxi- 
mation differs from approximations available at 
present. It differs from the 1/A^ expansion, since it 
keeps all the terms in N, from the Baym-Grinstein 
approximation as it still is of the <i>-derivable 
nature, from gapless approximations |20l l32l|. etc. 
Of course, at the Hartree-Fock level there is no 



"collisional" dissipation in the scheme (the Lan- 
dau damping may exist H^)- Nevertheless, the 
"collisional" dissipation can be introduced phe- 
nomenologically at the quasiparticle level, i.e. by 
adding a collision term to the Vlasov kinetic equa- 
tion resulted from this modified Hartree-Fock ap- 
proximation. Such an approach looks preferable, 
e.g., for the DCC simulations as compared 
to that based on the Baym-Grinstein approxima- 
tion. The "collisional" dissipation naturally ap- 
pears in higher-order <I>-derivable approximations, 
e.g. through the sunset self-energy diagrams. The 
corresponding modifications of these higher-order 
approximations in the way described here are much 
more involved and will be reported elsewhere. 
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Appendices 



APPENDIX A: TADPOLE 



In this appendix we present the limiting cases of the tad-pole integral introduced in Eq. H3.20|l : 

where e = (fc^ + Apy/'^/T. For T <^ M 

Q^iM, T ^0) = (27r)-3/2T3/2Mi/2 cxp(-M/r), (A.2) 



while for M <^T 



\M ^ 0,T) = 



1 1 M 1 M 
12 ~ 4^Y ^ 8^ V T^j 



(A.3) 



APPENDIX B: SECOND DERIVATIVE OF CJT EFFECTIVE POTENTIAL 

The equations of motion (|3.16() - (|3.18() in fact are conditions which determine the extremum of the CJT 
effective potential. However, stable solutions correspond only to minima of this potential. To distinguish 
between minimum and maximum of the potential, the sign of the second derivative of the potential over 0, 
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under the condition that the Green functions solve the gap equations with this (j), should be determined. 
The first derivative leads to 

dFgHF(</',G(^)) _ dV^M<P^G) , 5l^gHF(0,G) dG _ . 



<I^Mi~H, (B.l) 



d(j) d<j) 5G dq 

since 

oG 

due to equations of motion for Green functions. To obtain Eq. IjB.ip we have used the equations of 
motion (|XT7|) and (|XTB|) . Note that dVgHF{(t},G{(j)))/d(f) = 0, if is a solution to Eq. (jXTB)) . Taking the 
next derivative one arrives at 

dV,HP(0GW)^ .^^^^ (B.3) 
dtp'^ d(p 

In order to find dM'^/d(j) we can use the equations of motion (|3.16|) and 13.17|l . Taking derivatives over 
(j> in Eqs H3.16|l and (|3.17|) and then resolving them with respect to dAl'^/dcf), we arrive at 

dMl ^ [M, + ^2{N + 2)R,] 

# [A'U l)i?.] + ^Ra [m - 5)AU - ^2{N -1)(N + 2)R^] 

(B.4) 

with 

To avoid uncertainties in the calculation of dQ^ {M,T) / dM at M <^ T, it is reasonable to directly use 
Eq. ifOll . ^ 
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